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Abstract A monomer-dimer model with a short-range attractive interaction 
favoring colinear dimers is considered on the lattice Z 2 . Although our choice 
of the chemical potentials results in more horizontal than vertical dimers, the 
horizontal dimers have no long-range translational order - in agreement with 
the Heilmann-Lieb conjecture [10]. 
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Introduction 

A liquid crystal, at low temperatures, displays a long-range order in the 
orientation of its molecules, while there is no complete ordering in their posi¬ 
tions. In this paper we present a model characterized by these two features. 
In particular we consider a monomer-dimer model on the two-dimensional 
lattice Z 2 characterized by different chemical potentials for horizontal and 
vertical dimers (/q, > q v to fix ideas) and by a short-range potential J > 0 
that favors collinear dimers. We prove that when the parameters satisfy 

5 

Hh > — J and n v < - -J, (0.1) 

the system has the properties of a liquid crystal. 

Onsager [14] was the first to propose hard-rods models in order to explain 
the existence of liquid crystals. In 1970 Heilmann and Lieb [8,9] studied sys¬ 
tems of monomer and dimers (hard-rods of length 2) interacting only via the 
hard-core potential, and proved the absence of phase transitions in great gen¬ 
erality. Then in 1972 they [10] proposed two monomer-dimer models (named 
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I and II) on the lattice Z 2 , where short-range attractive interactions among 
parallel dimers are considered beyond the hard-core interaction. Heilmann 
and Lieb claimed that these systems are liquid crystals. In particular they 
proved the presence of a phase transition, by means of a reflection positivity 
argument: at low temperature there is orientational order. Moreover they 
conjectured the absence of complete translational ordering for their models. 
A proof of this conjecture for the model I was announced in [10] by Heil¬ 
mann and Kjaer, but never appeared. Letawe, in her thesis [12], claimed to 
prove the conjecture by cluster expansion methods, even if the result has 
never been published in a journal. Letawe’s polymers are built starting from 
contours and the major difficulty seems to arise when she has to deal with 
a polymer lying in the interior of another one: the two polymers would not 
be independent. To overcome this problem, ratios of partition functions with 
different (horizontal or vertical) boundary conditions Z v /Z h are introduced, 
but it is not proved that these ratios are sufficiently small to guarantee the 
convergence of the cluster expansion. 

Numerical simulations related to the Heilmann-Lieb conjecture are per¬ 
formed in [15]. We also mention that, in absence of attractive interaction, 
systems of sufficiently long hard-rods were proved to display a phase transi¬ 
tion and behave like liquid crystals by Disertori and Giuliani [3], using a two 
scales cluster expansion and the Pirogov-Sinai theory. 

In the present paper we study a model obtained from the model I of 
Heilmann and Lieb [10], but while they suppose 

Mh = Hv =: fJ. and fi > — J , (0-2) 

we assume very different horizontal and vertical potentials as in (0.1). This 
choice of the parameters allows us to work with cluster expansion methods, 
by defining our polymers starting from regions of vertical dimers, instead of 
contours. The cluster expansion method permits to rewrite the logarithm of 
the partition function of a polymer system as a power series of the polymer 
activities. This expansion entails analyticity results and simplifies consider¬ 
ably the study of the correlation functions, which can be expressed in terms 
ratios of partition functions. Clearly the cluster expansion cannot hold in gen¬ 
eral on the whole space of parameters: it converges only when the polymer 
activities are small enough to compete with the entropy. A rigorous study of 
the conditions of convergence dates back to [6,7,16], by means of Kirkwood- 
Salsburg type of equations. In this paper we use a criterion proposed by 
Kotecky and Preiss [11] in 1986. Afterwards this criterion was compared to 
the previous ones, was improved and simplified in [1,2,4,5,13,17] (for a clear 
and modern treatment we suggest for example the last work). 

The paper is organized as follows. In the section 1 we introduce the model 
and we state the main results about its liquid crystal properties. In the 
section 2 we show how to rewrite the partition function as a suitable polymer 
partition function, following in part the ideas of [12]: our polymers turn out 
to be connected families of regions of vertical dimers and lines of horizontal 
dimers and monomers. In the section 3 we prove that the Kotecky-Preiss 
condition for the convergence of the cluster expansion is verified when the 
parameters satisfy (0.1) and the temperature is sufficiently low. Finally in 
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the section 4 we use the previous sections to prove the results stated in the 
section 1. The appendix A contains the study of a 1-dimensional monomer- 
dimer model, that is needed in the section 2. For the sake of completeness, 
in the appendix B we state the general results of cluster expansion needed 
in the paper. 


1 Definitions and Main Results: the Model and its Liquid Crystal 
Properties 

A monomer-dimer configuration on Z 2 can be represented by a bonds 1 
occupation vector a £ {0,1} B( A ) with hard-core interaction, namely: 

< 1 VxGZ 2 . (1.1) 

y~x 

If cc( Xj y) = 1, we say that there is a dimer on the bond (x,y) , or also that 
there is a dimer at the site x\ if instead 0 ( Xj y) = 0 for all y ~ x, we say 
that there is a monomer on the site x . Dimers on Z 2 may have two different 
orientations: vertical ( v-dimers ) or horizontal ( h-dimers ), according to the 
orientation of the occupied bond 2 . The model studied in the present paper 
favors one orientation of the dimers (the horizontal one), both via a chemical 
potential and via a short-range imitation. 

Let A be a finite sub-lattice of Z 2 . Consider a horizontal boundary condi¬ 
tion 3 , namely we assume that every site of Z 2 \ A has a h-dimers (with either 
free or fixed positions). Denote by the set of monomer-dimer configura¬ 
tions on A (we allow also dimers toward the exterior 4 ) which are compatible 
with the selected horizontal boundary condition. 

The Hamiltonian, or energy, of a monomer-dimer configuration is defined 
as 

._ J .#/sites of A with "I , /A — Mv »r sites of A with "l < 

A ' 2 l monomer J 2 tv-dimer / 

J / I" sites of A with h-dimerl ("sites of A with v-dimer'l \ 
+ — I if < but h-neighbor also to a > + ffi but v-neighbor also to a > ) . 

^ V Iv-dimer or a monomer J Ih-dimer or a monomer J J 

(1.2) 

We assume that the parameters appearing in the Hamiltonian satisfy 

Mh A J , ■ ' fiy , J > 0 . (1*3) 

1 Two sites x = (*h,a;v), y = (yh,y v ) £ Z 2 are neighbors (x ~ y) if \xh — j/h| + 

\x v — J/v| = 1 • A pair of sites ( x,y ) is a bond if x,y are neighbors. B(Z 2 ) denotes 
the set of bonds. 

2 Two sites x = (xh,x v ), y = [yh,y v ) £ Z 2 are h-neighbors if x v = y v and 

\xh — 3/h| = 1, they are v-neighbors if Xh = yh and |x v — y v \ = 1. A bond (x, y) £ 

B{1?) is horizontal if x, y are h-neighbors, it is vertical if x,y are v-neighbors. 

3 The external boundary of A is d ext A := {x £ 1? \ A \ x neighbor of y £ A}. The 
internal boundary of A is instead <9/1 = d mt /l := {x £ A \ x neighbor of y £ 1? \ A}. 
We set A := A U <9 ext /l. 

4 Namely we allow dimers having one endpoint in A and one in 1? \ A. 
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In this way, if the horizontal boundary condition with free positions is cho¬ 
sen 5 , then the ground states in (i.e. the configurations minimizing the 
energy under the given condition) are exactly the configurations where every 
site has a h-dimer. The partition function of the system is 

Z\ := ^ e - 0HA(a) (1.4) 

a£®!; 

where the parameter /? > 0 is the inverse temperature. 

Remark 1.1 We want to show that the Hamiltonian (1.2) essentially corre¬ 
sponds to the model I introduced by Heilmann and Lieb in [10], except for 
the important fact that we allow the horizontal and vertical dimer poten¬ 
tials fi h, n v to be different, while they take /Zh = /u v = T ■ We can introduce 
another Hamiltonian (that maybe is written in a more natural way; see fig.l): 

Hj i := — 4 % ^{h-dimers in A} — Mv ^{v-dimers in yl} + 

_ j gif pairs of neighboring\ (1-5) 

’ \ col i near dimers in A j 




Fig. 1: The same monomer-dimer configuration on the lattice A and the corre¬ 
sponding energies in accordance to the Hamiltonian (1.2) (on the left) and to the 
Hamiltonian (1.5) (on the right). A horizontal boundary condition is drawn in grey. 


The monomer-dimer model I in [10] is given by the Hamiltonian (1.5) with 
Th = Tv = T- when A is a rectangular lattice of even sides lengths with 

5 Also fixed positions work, provided that the positions of the two h-dimers at 
the endpoints of each horizontal line of A allow a pure dimer configuration on that 
line. 
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periodic boundary conditions (torus). It is easy to show that when ff is a 
torus the two Hamiltonians (1.2), (1.5) describe the same model; indeed they 
only differ by an additive constant which does not affect the Gibbs measure: 


H a + \A\ = H a 


( 1 . 6 ) 


since 


| - 2#{h-dimers in ff} = |ff| - with } = 

{: 


_ ///sites in ff with) 
" \ monomer / 


i al / sites in ff with \ 
' "tv-dimer / 


OM-S ]■ • it — 44 /sites in ff with/ . 

2#{v-dimersmff} - #( v _ dimer ) , 


, , ( sites in ff with h-dimer (v-dimer)) 

|ff | - 2 P a . irs of neighboring I ni # J and h _ neighbo r (v-neighbor) to 

v J (^another h-dimer (v-dimer) J 


__ ///sites in ff with/ 
" \ monomer / 


{ sites in ff with h-dimer (v-dimer) > 
and h-neighbor (v-neighbor) also 
to something different 


On the other hand when ff has horizontal boundary conditions the two Hamil¬ 
tonians (1.2), (1.5) are not exactly equivalent. Indeed it holds 6 


Ha + 


/t h + J 
2 


14 + 


— 44 /sites in <9/ nt ff\ _ jj 
2 ” /without h-dimer J 


(1.7) 


when the following conventions are adopted in the definition (1.5): if only 
half a dimer is in ff while the other half is in Z 2 \ ff, it counts if only one 
dimer of a pair of neighboring colinear dimers is in ff, while the other one is 
in Z 2 \ ff, this pair counts 


The monomer-dimer model that we have introduced, in a certain region 
of the parameters corresponding to large horizontal potential, small vertical 
potential and low temperature, behaves like a liquid crystal. This means that 
the model exhibits an order in the orientation of the molecules (dimers), 
while there is no complete order in their positions. 

The following results will give a precise mathematical meaning to these 
statements. First we introduce some observables attached to the sites, asking 
questions as “Is there a horizontal dimer at site xT ’, “If so, is it positioned 
to the left or to the right of aTo measure the absence or presence of 
some kind of order, at a microscopic level we study the expectations and 
the covariances of these quantities according to the Gibbs measure, while at 
a macroscopic level we introduce a suitable order parameter and study its 
expectation and possibly its variance'. 

6 d v , dh denote respectively the vertical , horizontal component of the bound¬ 
ary; e.g. <9 v ff := {x £ A | x h-neighbor of y £ Z 2 \ ff} and <9hff := {x £ 
A | x v-neighbor of y £ Z 2 \ ff}. 

' When the expectation of the order parameter is zero but the variance is not, a 
small perturbation can lead to a spontaneous order of the system. 
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Define the following local observables 8 

fh,x ■= l(a; has a h-dimer) , f ViX := l(x has a v-dimer) ; (1.8) 

f\,x '■= l( aT has a left-dimer) , f rx := \[x has a right-dimer) . (1.9) 

Clearly f^ x = f\ x + f rx and x + f vx < 1. In the following we denote the 
Gibbs expectation of any observable / by 

(f)\ ■= Jh E /(«)e-^ (a) . 

A ae 

We denote by N the minimal distance between any two vertical compo¬ 
nents of the boundary of A and our only assumption on the shape of A is 
that N — i oo as A /* Z 2 . To fix ideas one could think that A is a rectangle 
(in this case N would be simply its horizontal side length), but actually we 
will need to consider also non-simply connected regions. 

There exists fio > 0 depending on //h, Abo J only and No(fi) depending 
on /3, Hh, J only such that the following results hold true. 


Theorem 1.2 (Microscopic expectations) Assume that J > 0, /./h + J > 
0 and 2p v + 5J < 0. Let (3 > (3 q. Let A C Z 2 finite having N > Nq(/ 3). Let 
x £ A such that disth(a;, dA) > N 0 ((3). Then 


<Ax>n > \ 


(1.10) 

As a consequence: 


(A ,x)\ > 

(1.11) 


\{fi,x)\-(hx)\ \ < 2 c . 

(1.12) 


Theorem 1.3 (Microscopic covariances) Assume that J > 0, /Jh + J > 0 

and 2 fi v + 5J <0. Let /3 > 0o- Let A C Z 2 finite such that N > Nq((3). 
Let x,y G A such that disth(a;, dA) > N 0 (/3), disth(y, dA) > N 0 (fi) and 
disth {x,y) > Nq(/3). Then: 

^ e-T(«^(*.»)-i) , (1.13) 

|</r,x fr,y)\ {fr,x)\ ^ (*,,)-!) > (1 . 14) 

\{h,x fr,y) h A - {fl,x) h A {fr.y)'\' < ^ e“^(,,y)-l) . (1 . 15) 

The definition of m is clarified in the Appendix (lemma A.5); anyway it can 

Q / / \ \ 

be sufficient to know that m = e~“ 2 (1 -)- o(l)) as fS —► 00 . 

8 We say that the site x has a left-dimer if there is a dimer on the bond (x, x — 

(1, 0)) , a right-dimer if there is a dimer on the bond (x, x + (1, 0)) . 
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The density of lattice sites occupied by h-dimers/v-dimers is respectively: 

^ := uf E A,® ’ := \X\^ A,®' (i- 16 ) 

' ' xGA ' ' xeA 

A parameter measuring the orientational order of the dimers is 

^orient. * = ■ (1*17) 

Corollary 1.4 (Orientational Order Parameter) Assume that J > 0, 
Hh + J > 0 and 2/z v + 5J <0. Let ft > /3q. Let Ac Z 2 finite, having 
N > 2N 0 (/3) . Then 

(A mient .)\ > (l - 2 M) (1 - 4c . (1.18) 

Hence 

lim lim inf (Adent.)* = 1 • (1-19) 

P/'oo A /*Z 2 


The corollary 1.4 shows that fixing ft sufficiently large and then choosing A 
sufficiently big (more precisely the distance N between vertical components 
of dA must be large enough), the average density of sites occupied by h- 
dirners is arbitrarily close to 1: in other terms the system is oriented along 
the horizontal direction. 

The majority of sites is occupied by h-dimers. But there can still be some 
freedom, indeed we may distinguish the h-dimers in two classes according to 
their positions: a h-dimer is called even (resp. odd) if its left endpoint has 
even (resp. odd) horizontal coordinate. The density of lattice sites occupied 
by even/oclcl h-dimers is respectively: 


1 _ 2 _, 

— \(x has an even h-dimer) = — f I X , 

I I xeA I I „<=/! 


xu even 


j _ 2 _^ 

J'odd : = p| l(z has an odd h-dimer) = yj\ A* • 


xeA 


A 


xu even 


A parameter measuring the translational order of the h-dimers is 

"Atransl. •— ^even ^odd • 


( 1 . 20 ) 


( 1 . 21 ) 


Corollary 1.5 (Translational Order Parameter. Part I) Assume that 
J >0, Hh + J >0 and 2/r v + 5 J < 0. Let ft > fio- Let A C Z 2 finite such 
that N > 2N 0 ((t). Then 


(Aransl .)\ < (1 — 2 


Noifi) 

N 


2 e“ 


, r, N o(P) 
+ 2 A\A 


( 1 . 22 ) 


Hence 


lim lim sup I (Aransl.) \ 


= 0 . 


(1.23) 



Corollary 1.6 (Translational Order Parameter. Part II) Assume that 
J >0, /ith + J >0 and 2/z v + 5J < 0. Let f3 > (3 q. Let A c Z 2 finite such 
that N > 2 7Vo(/3) . Then 


<(* 


transl. 



((Aransl.)A ) 2 < 


1 9 TO 

R (l-e"T)2 


W) 

N 


Hence for fixed (3 > f3 d 


(1.24) 


lim ((Z\ 

yl/'Z 2 V 


transl. 



«AranBl.)J0 2 = 0 . 


(1.25) 


The corollaries 1.5, 1.6 show that fixing (3 sufficiently large and then 
choosing A sufficiently big (in particular the distance between different com¬ 
ponents of d w A must be big enough), the mean value and the variance of 
the difference between the density of even h-dimers and the density of odd 
h-dimers are arbitrarily close to zero. In other terms, at large but finite /?, 
there is not a spontaneous translational order for the h-dimers. 


Remark 1.7 The bounds (1.22) hold for any kind of horizontal boundary 
conditions, but in some particular cases it is possible to obtain a better 
result by a symmetry argument. Assume that A is a rectangle with N + 1 
sites in each horizontal side. If N + 1 is odd, by choosing horizontal dimers 
with free positions at the boundary one obtains 

(AranslR = (z'everR - (^odcR = 0 (1.26) 

for all parameters (3, J, //>,, p v . To prove it consider the reflection on A with 
respect to the vertical axis at distance ^ from d w A: this transformation 
induces a bijection T : It is easy to check that Ha(T(ol)) = 

Ha{u) , i4ven (T(a)) = n odd (a ), v 0 dd{T(a)) = i/ eve n{a) for all 
On the other hand if N + 1 is even, by choosing periodic boundary conditions 
one still obtains 

(Aransl.)r' = 0 (1.27) 

for all parameters (3, J, /.%, /J v ■ To prove it one can consider the reflection 
on A with respect to two vertical axis at distance from each other: it 
induces a bijection from S>^ el ' to itself having all the previous properties. 


2 Polymer Representation 

In this section we show how to rewrite the partition function Z\ as a 
polymer partition function of type (B.l). This representation will be suitable 
for applying the cluster expansion machinery (see Appendix B) in a regime 
of large horizontal potential, small vertical potential and low temperature. 

We start by isolating the “few” vertical dimers. Associate to each monomer- 
dimer configuration a £ the set 

V = V (a) \= {x £ A \x has a v-dimer according to a} . 
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Partition V into its connected components (as a sub-graph of the lattice 9 
Z 2 ): 

n 

V = U Si ’ y A ^i , dist Z 2 (Si, Sj) > 1 Vi ^ j 

i=i 

where the family SZ A is defined by 

S € y A ^ S C A , S connected (as a sub-graph of Z 2 ) , 

every maximal vertical segment of S has an even number 
of sites , 

S does not contains those sites of d^A that necessarily 
have a h-dimer because of the boundary conditions. 

( 2 . 1 ) 


The knowledge of the set V (or equivalently of Si ,..., S„) does not determine 
completely the configuration a of the system, since on A \ V there can be 
both h-dimers and monomers. Anyway a fundamental feature of the model is 
that the system ond\P can be partitioned into independent 1-dimensional 
systems. Introduce the family JZ A (V) defined by 


L € Jz? A (V) L is a maximal horizontal line of A \ V . (2.2) 


The Hamiltonian (1.2) rewrites as 


Ha = E |Si| + \ |9h5<l 

/ih + J „ f sites of L 
2 " \ monomer 



- I d v Si n 8A\ 


withl 


J 

2 


# 


{ 


sites of L with h-dimer 
but h-neighbor also to a 
monomer or to Ui Si 



Hence the partition function (1.4) rewrites as (see fig.2) 


7 h 


= ^ ^ e -p( tlL 2 ty -\Si\ + i\d h s i \ + ±\dvS i ndA\') j-j- 

Le&A(ViSi) 


n> 0 


s u ...,s n e^A t=i 

dist(Si,Sj)>l 


( 2 ; 3 ) 

where Zl is the monomer-dimer partition function of the line L, consid¬ 
ered as a sub-lattice of the 1-dimensional lattice Z, with suitable boundary 
conditions: 


Zl := E g &Hl (cX-l) gl\,x\ (c*a:i ) g-^r ,x r ( a x r ) ^ 

OtL^^L 


9 On any graph the distance between two objects is defined as the 
length of the shortest path connecting them. In particular dist Z 2 (5, S') := 

infugs.^es' dist Z 2 (x,y) for all S,S' C Z 2 and dist Z 2 (*, y) := \xh - t/h| + \x v - y v \ 
for all x = (x h ,x v ), V = (Vh, y v ) € Z 2 . 
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Fig. 2: A monomer-dimer configuration on A and the corresponding regions 
S'i,S , 2 ,S , 3 and lines Li,...,Lis £ JzfA(UiSi). Given the positions of the regions, 
the configurations on the lines are mutually independent: the arrows represent the 
energy contributions of type J/2. A horizontal boundary condition is drawn. 


An explanation of the notations introduced in (2.4) is required. £>l denotes 
the set of monomer-dimer configurations on L (dimers can only be horizontal, 
external dimers at the endpoints of L are allowed); 


H l 


Fh + J 
2 


*{ 


sites of L with 
monomer 


+ 



sites of L with dimer) 
but h-neighbor also to ; 
a monomer J 


xi, x r denote respectively the left, right endpoint of the line L (which even¬ 
tually may coincide): observe 10 that because of (2.2) 


(J xfiL) = u (WlXUiW), (2.5) 

L£se A (UiSi) ' ' ' ' 

U xfiL) = U (2.6) 


L£^ A (UiSi 


10 di, d r denote respectively the left, right component of the vertical boundary; 
e.g. d\A := {a: £ A \ x — (1,0) G I? \ A} and d v A := {x £ A \ x + (1,0) £ I? \ A}. 
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finally 11 

if x\ G U,; d^Si => 7i iXl := (-00 -/if 0) 

if xi G d\A, on X\ — (1,0) it is fixed a 1-dimer =t- I\ tXl := (—00 0 —/3f) 

if xi G d\A, on xi —(1,0) it is fixed a r-dimer =» I\^ XI := (0 —00 — 00 ) 

if xi G d\A , on xi —(1,0) there is a free h-dimer => I\_ Xl := (0 0 — /if) 

(2.7) 


and, similarly, 

if x r G Ui 9f xt 5j =t- I I>Xl := (-/if -00 0) 

if x r G d x A , on x r + (l, 0) it is fixed a r-dimer => 7 r)Xt := (0 —00 — /if) 

if x r G 9 r yl, on x r + (l, 0) it is fixed a 1-dimer => I I Xr := (—00 0 — 00 ) 

if x r G d r A , on x r + (l, 0) there is a free h-dimer => / r>a , r := (0 0 —/if) • 

( 2 . 8 ) 

The 1-dimensional systems described by Zl, L G ,2a (UjS)), are studied in 
the Appendix A. 

In the form (2.3) of Z \, the weight of the regions (Si ,..., S n ) is not a 
product of the weights of each region Si, because of the lines L connecting 
different regions. Therefore the regions Si G ^2a are not a good choice for a 
polymer representation of the model. In order to decouple some regions from 
some other ones, it is possible to do a simple trick. It is convenient to deal 
in different ways with the endpoints lying on d ext Si and those on dA; hence 
given a line L G Jz^UiS)) we set 

£i,*i := l(xi G (Uj 5® xt Sj) D A) , r/i^j := 1- £i, Xl ( = 5) l(xi G (<9i A) \ U* d\Si) 

£i, Xl ■= 1 (^r G (Uj <9f xt 5i) (~l A) , r/ TtXi := 1- £ I:Xi (2 = 6) l(x r G (d r A) \ Uj dr Si) 

Using the notations of the Appendix A, given a line L G .^(UiS)) we in¬ 
troduce the two vectors representing the boundary conditions outside its 
endpoints X\, x r : 



/ e lT ’ x A l ) \ 

B\,xi '■= e /3 ** h 4 + A*i( m )^ 

, B,, Xi := e^Ar) 

\ e -/3 a ^ + V,x r ( m ) / 


then to shorten the notation we set 
1 


^l,xi * — 


Tv Bui E '" 


TV 


e[ x) b^ Xi 


11 The possible states of a site x G L are three: “1” =left-dimer namely a 
dimer on the bond (x,x — (1,0)), “r v =right-dimer namely a dimer on the 
bond (x,x + (1,0)) , “m”=monomer. Here we think I\ tXl , I r , Xr as vectors: I\ )X , = 
(7],xi(0 I),xi (r) 7i,x, (rn)) and 7r,x r — (7r,x r (0 A,x r (0 7r,x r (m)). 
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Now define 


Rl 


Zl 


A^ 1 b m ’ 




- c * 


i^i 


(2.9) 


and, using 2z? as an abbreviation for S£a (U,;5j;), rewrite the quantity YIlg& 
by means of elementary algebraic tricks: 



By identities (2.5), (2.6) it holds 


II Cfe- = 


Leif 


n 

Leif\jr 


n ^ 
n 


^xG(Ui^ xt 5i)n2l 
supp Jff 



n 


x€d T A\Uid r Si 


n 


xG(Ui^f xt 5i)nyi 
x(£ supp 


By substituting into the previous formula and thinking J(T = {L\,... ,L p }, 
we find out 12 


n 

Leif 




■E = £ 

p>0 ^ Li,...,L p eif 
L h ^L k Vh^k 



n '' 


* ®e(Ui9; xt Si)na 
x^UkLk 


( 2 . 10 ) 


12 In the first product on the r.h.s. of (2.10) the shorten notation b\/ TtX means: 
take b\, x if x £ d\A, take b T , x if x £ d r A; notice that d\A and d r A are disjoint for 
N > 1. In the last product instead the shorten notation b T /i tX means: take b 1}X if 
x £ df xt Si only, take bi, x if a: £ d^ xt Si only, and take the product b r , x b\ tX in the 
case that x belongs to both d\ xt Si and d x xt Sj. 
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Now substitute (2.10) into (2.3), using also the fact that |/1| = ^” =1 |5,| + 
Eie^OJiSi) l L l’ and obtain: 



■Es e n 

n>0 5 i ,...,S n 6 *=1 

dist (Si, 5j) > 1 Vi^j 




js.l 


n 


xGd v ADd v Si 



•e^ e fn^) ( n m.*) 

P> o Li,...,ip€^U(UiSi) \fc=l / V a:e(Ui9; xt Si)nyl / 
Lk^LhVk^h x^UkLk 


( 2 . 11 ) 


The next step is to partition |J " =1 Sj U ljfc=i Lk into connected components 
as a sub-graph of Z 2 , where Z 2 is the lattice obtained from Z 2 by removing 
all the vertical bonds incident to the lines L k : 


n p q 

U^U Lk = U SU PP p t 

i=l k =1 t=1 


Pt G ^ Vi , distg 2 (suppP t ,suppP s ) >1 Vi 7 ^ s 

where the family (yes, it is finally our family of polymers! see fig.3) is 
defined by: 


& A := 


{P = (( Si)i£i , (Pfe)feeic) 
(5j)iej G ^ I 


(5i)i G (L fc )fc e ^(UiSj)} , 

( 2 . 12 ) 

0 < |/| < oo 

Si G &a V* (2.13) 

distz 2 (S,, Sj) > 1 V* ^ j , 


( Lk)keK G ^^(Uig/Sj) 


' 0 < \K\ < oo, \I\ + \K\ > 1 
L k G ^(UiSi) Vfc 
L k 7 ^ LhVk jtz h 

(U iSi) U (Ufcifc) connected in Z 2 . 

(2.14) 


The identity (2.11) now rewrites as 

Z A = C A J2 E f[sA(Pt) n^,P.) ( 2 - 15 ) 

q >0 Pi,...,P q e&A t-1 t<s 

by setting, for all P,P' G with P = ((Sj) ie j, ( L k ) keK ), 

C'a := A ^ 1 H 6 1/r) 


(2.16) 
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Si 
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Li 



















Fig. 3: The first three pictures represent three different examples of polymers P € 
S?a- The set represented in the last picture is not a unique polymer since it is not 
connected in Z 2 (even if it is connected ini?). 


Qa{P) 



/ e -fi( tSL 2^\Si\ + i\d h S i \) 


A 


I Si | 

1 




x £ U keK L k 


(2.17) 


S(P.P') 


1, if distg 2 (P, P') > 1 
0 , otherwise 


(2.18) 


The identity (2.15) finally shows that the partition function Z\ , up to a 
factor Ca, admits a polymer representation of the form (B.l). 

It is convenient to bound the polymer activity qa by a simpler quantity. 
Using the proposition A.9 plus the lemmas A.6, A.8 and the fact that |(9h<5.;| > 
2 , one finds: 


Qa{P ) < e(P) := 



(2.19) 


with the 7 l’s defined by the equation (A.9). 
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3 Convergence of the Cluster Expansion 

In the previous section we rewrote our partition function as a polymer 
partition function up to a factor Ca (see formula (2.15)). In this section we 
will find a region of the parameters space ph, /x v , J where the condition (B.2) 
is verified by our model at low temperature, so that the general theorem B.l 
about the convergence of the cluster expansion will apply to our case. 

Theorem 3.1 Assume that J > 0, /Zh + J >0 and 2/r v + 5J < 0. By 
choosing 

TY1 

a(P) := — | suppP| VPe^yi (3.1) 

the conditions 

E 9 {P)e a(P) <^ Vied, (3.2) 

Pe&A 

supp P3x 

J2 q{P) e a(P) < a(P*) VFe &> A (3.3) 

Pe&A 
5(P,P *)=o 


hold true, provided that /3 > /3o and N > No(/3) (N is the minimum distance 
between two vertical components of dA). Here (3 q > 0 depends on ph,/i v ,d 
only, while Nq(/ 3) depends on /3,/th, J only; they do not depend on A, P *, x. 

Corollary 3.2 Assume that J > 0, /% + J > 0 and 2^ v + 5 J < 0 . Suppose 
also that (3 > 0o and N > Nq(/ 3). Denote by ^SZa the set of clusters 13 
composed by polymers of SZ A ■ Then the partition function (1.4) rewrites as 

Z\ = Ca exp ( E* U 4( p t)t)) (3-4) 

\p t )teV3»A ' 

where we denote E*p t ) t e^ : = E g >o E (P t )? =1 e*0»A and 

Q 

U A (Pi,...,P q ) := u(Pi,..., P q ) n OA(Pt) ■ (3.5) 

t —l 

Remind that Ca is defined by (2.16), qa is defined by (2.17) and u is defined 
by (B.4), (2.18). Furthermore for all $ C SZ A it holds 

E* lM( p *)t)| < E l^( P )l ea(P) (3-6) 

(Pt)te^A Pes° A 

3t-.p t es Pe£ 


where a is defined by (3.1). 

13 As explained in the Appendix B, using the definition (2.18) for 5, a family of 
polymers (Pi ,..., P q ) is a cluster iff U® =1 supp Pt is connected in 1?. 
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Proof The corollary follows from the general theory of cluster expansion (the¬ 
orem B.l), since Z\ admits a polymer representation (2.15) and satisfies the 
Kotecky-Preiss condition ((3.3), \qa\<q)- □ 

For ease of reading, in the following of this section we will denote 

£* == £ s £ “ d £’==£* £ 

(Si)i n (Si)? =1 e&S"A (L k )k P y (L k )l =1 e^ A (UiSi) 


where hP^A^iSf) are the projections of the polymer set Pa defined 

in (2.13), (2.14). The next lemmas provide the entropy estimates that will 
be needed in the proof of theorem 3.1. 

Lemma 3.3 If UjS, ^ 0, namely n> 1 , then 

1 < 4£‘l s ‘l . (3.7) 

{L k )k 


Proof Fixp > 0 and denote by ^Jzf^(U,;£)) the set of (Lfc)fc=i € ^ 2 L5f y i(U,;S' i ). 
Given (Lfc)fc=i G (UiSi), each line L /. has at least one endpoint on 

U id% xt Si , since (UjSj) U (U kLk) have to be connected in Z 2 . Therefore the 
number of ways to choose each L is at most J2i |c?® xt ^ | < 2 |<?»|. Since 

the Lfc, k = 1,... ,p, must be all distinct, it follows that 


|^?j p) (U < S i )| < (2^. 15,1) (2^15,1-1) ••• (2^15,1-p+l) . 
Therefore 


£‘ 1 = £sl^“< u - s ‘>b£ 

{L k ) k V P 


^ Tli I'S'i 
P 


= 2 2 £i 


□ 


Lemma 3.4 Let x £ Z 2 . For all s > 2 

16 

#{5 C Z 2 connected | |Sj = s, S 9 x} < — 4 4s . (3.8) 

O 

Proof Given a connected graph G and one of its vertices x, there exists a walk 
in G that starts from x and crosses each edge exactly twice 14 . Therefore 


ff{S C Z 2 connected | |S| = s , S 9 xj < 

2s 

< Y2 #{S connected sub-graph of Z 2 | |edges of Sj = e , 5 9 i} 


e—s— 1 
2s 


< Y2 ff ^ walks in Z 2 that start from x and have lenght 2e} 

^4s+2 


e=s— 1 
2s 


< 4 2e < 


e—s— 1 


□ 


14 


This can be easily proven by induction on the number of edges. 
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Lemma 3.5 Let A C Z 2 finite. For all s > 2, 1 < d < oo 

32 

#{S C Z 2 connected | |S| = s, disth(<5, A) = d} < — |A|4 4s . (3.9) 

O 

Here disth(>S', A) := inf^gs^gA disth(a^,j/) and the horizontal distance between 
x = (sh,a;v), y = ( yh,Vv ) G Z 2 is defined as 


disth (x,y) 


\xh~Vh\ if x v = y v 
+oo if x v y v 


(3.10) 


Proof Observe that disth {S, A) = d if and only if there exists a horizontal 
line L, \L\ = d + 1, having one endpoint on d w A and the other one on d v S. 
Therefore: 

#{S C Z 2 connected | 15*1 = s, disth(<5, A) = d} < 

< 'y ' ff{S C Z 2 connected | 15 1 1 = S , d v S 9 other endpt. of + } 

L horiz. line, |L|=d+l, 
d w A3 one endpt. of L 

< 2\d w A\ #{S C Z 2 connected | |5| = s, S 9 0} < 2\A\ ^ 4 4s . 

O 


For the last inequality we have used the lemma 3.4. □ 

Lemma 3.6 Let n > 1 . Let T be a tree over the vertices {1, ... , rz} . Let 
Si > 2 for all i = 1,..., n and dij > 2 for all (z, j) € T■ 

Then given A C Z 2 and 1 < d < oo 


#{($)?= i e | dist h (Si, A) =d, | Si | = Si Vz, 

dist h (S), Sj) = dij\/(i,j)&T} < 


i =1 


< Ain(f 4l '‘*-' ,T( " 


while given x £ Z 2 

#{(50?=1 G I Si 9 X, I Si I = Si Vz 


disth(5i,5,-) = dyV(z,j)GT} < 


< n(f 4 4s ‘ S i esrW 


(3.11) 


(3.12) 


Here deg-^z) denotes the degree of the vertex i in the tree T■ 

Proof Let start by proving the inequality (3.11) by induction on n. If n = 1, 
then the tree T is trivial and (3.11) is already provided by the lemma 3.5. 
Now let n > 2, assume that (3.11) holds for at most n—1 vertices and prove 
it for n . It is convenient to think that the tree T is rooted at the vertex 1 and 
denote by z the relation “vertex j is son of vertex i in 7”” and by T (z) the 
sub-tree of T induced by the vertex i together with its descendants. Then, 
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denoting by -/Vr,i (^i ( s i)ieT) (^ii)(i,j)er) the cardinality on the l.h.s. of 
(3.11), it holds 


Nt ,i (A d; (si) ie r, (dij)(ij)eT ) — 


- E 

Sie^A, |Si|=ai 
disthC-Si ,A)=d 


n n t(v),v(si, dlv'i ( s i)ieT(w)) {dij)(i,j)£T{v)) ■ 


Since T(v) has at most n — 1 vertices, the induction hypothesis gives 


NT(v),v(Sudlv',(Si)i£T(v)>(dij)(i,j)£T(v)) < S 1 JJ fy ^ s i Sr(v,< ' ^ . 

ieT(v) ' ' 


Then by substituting in the previous identity, bounding deg-^,,) (i) by deg 7 -(z) 
and using the lemma 3.5, one obtains: 


-^T,l d’, , {dij E |^1| 11 

zer 


— 4 4s i „ de Sr(®) ) 

3 i ) 


This concludes the proof of (3.11). 

In order to prove the inequality (3.12), denote by N (- : (cc; (sj)j 6 7 -, (dii)(ij)er) 
the cardinality on the l.h.s. of (3.12) and observe that 


d^T, l( a '» ( s i)ieT? ( dij)(i t j)£T ) 


E 

Sie^A, |5i|=si 
S i 3 x 


LI NtmASu div’i ( s i)ieT(v)i (dij)(i,j)eT(v)) ■ 

V4—1 


By (3.11) we already know that 


N r(v),v(Si,di v ;(si) ie T(v),(dij)(i,j)eT(v)) E «i ( y 4 4 s ‘s i St(,,,( ) j . 


i£T(v 


Then by substituting in the previous identity, bounding deg-y-^ (i) by deg-^z) 
and using the lemma 3.4, one obtains: 


d^T, l(®> ) (dij ) (i.j 


1 er) E 


n 

ter 


— 4 4s i 

3 


deg r (i 


which proves (3.12). □ 

Proof (of the theorem 3.1) According to the definition (2.18), the condition 
S(P,P*) = 0 implies that supp P D [suppP*]i ^ 0, where [A]i := {x £ 
Z 2 | dist^a ( x , A) < 1} . Therefore 


£ Q(P)e a(P) 

pe&A 

6 (P,P* )=0 


< £ £ Q(P)e a(P) 

ai£[supp P*]i pz&a 

supp P 3 x 


< 4 I supp P* I max 

x£A 


£ Q(P)e a{P) - 

P££? a 

supp P 3 x 
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Thus, by choosing a(P) := "f | supp P\ for all P G LPa, the inequality (3.3) 
will be a consequence of (3.2). 

We have to prove the inequality (3.2). It is worth to write down explicitly 
the quantity we will work with (see the definitions (2.19) and (3.1)): 


g(P)e a{p) = 



IT 

k =1 


> 2 I ^ k I 


7 L k 


for all P G LPa, P = ((<Si)" = i, (Tfc)^ =1 ) • Notice that if suppP 9 x, the site 
x may belong either to a region S, or to a line Lk ; hence we can split the 
sum on the l.h.s. of (3.2) into two parts: 


Y e(P) e a(P) = Si + £2 (3.13) 

supp 


with 


* - E‘ 


) y ri e 

(Si)i 

UiSi3x 

Vi 


a - E‘ ( 


Y 

(Si)i ' 

V i / 

( Lk)k k 


(3.14) 


(3.15) 


During all the proof o(l) will denote any function u> = u>(/3, /i h, J) such that 
u> —> 0 as /3 — > 00 and u> depends only on (3, /./>, , J (in particular it does not 
depend on the choices of A C Z 2 , x € Z 2 , P G &a). 


I. Study of the term S\. 

We hx a family of regions (.S'j )" =1 that contains the point x; we also 
assume that ,^_5?/i(U.;>S',;) is non-empty, otherwise the contribution to is 
zero. By the lemma 3.3 it holds 

e ~it\ L k\ry Lk < 4^il s d max JJ e - ^ 1 ^ 7 i, fe (3.16) 

(Lk)k 



where the maximum is taken over all (Lk)k G ^L5fyi(Uj5j). The factor 7 ^ 
can take two values (see formula (A.9)), both smaller than 1 for /3 sufhciently 
large (uniformly with respect to Lk), since each line Lk must have at least 
one endpoint on to ensure that (U»Sj) U (Lh-Lfc) is connected in Z 2 . 

Obviously n > 1 in order for U” = 1 >f>j to contain the point x. It is convenient 
to consider separately the case n = 1 and the case n > 2 : 

Si = S[ + S'{ . 
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The case n = 1 is easy to deal with, simply by bounding the r.h.s. of (3.16) 
by 4J S I and using the lemma 3.4. Precisely: 


^1:= E 


r) \S\~PJ 


e ” |ifel 7 L k 


s&y A 

S3x 




< 


E 


se^A 

S3x 


“ ^ 3 

s> 2 

16 


?) 1^1 ^|*S'| 

_ tt ) S ~PJ ^s 


(3.17) 


= ~tt 4 10 e _/3 (Mh_Aiv+J) (1 + o(l)) . 

O 

Now assume n > 2. Fix a family of lines (T/c)fe=i € ^ 2 LSf / i(Ui5',). We 
can consider the graph G = G((Si)i, {L k )k) with vertices * £ {1,..., n} and 
edges k £ {1,... ,p} : the edge k joins the two vertices *, j iff the line L k has 
one endpoint on <9® xt .S) and the other one on d° xt Sj. In the graph G there 
can be multiple edges, loops and pseudo-edges with a single endpoint. The 
graph G is connected (it follows from definition 2.14), hence G admits at 
least one spanning sub-tree T. And clearly, since each factor e~^ Lk ^ r )L k is 
smaller than 1, 


n 

k=l 


e ™ l ifc l < 


n 

feer 


e " 2 ' |Lfc| 7 L k < 


n 

(ij)er 


e -^(dist h (5i,Sj)-l) s 


where 7s,S' := (|e + e ( dlsth (S,S ) i)^ ^ _ Therefore: 


max 1 \e 2,ljfe| 7 l, < max 
(L k ) k 11 


,-f lifcl 


n 


— ^(disth(S»,Sj)-l) 


T tree over 

(ij)eT 


e 2 


75i.fi,- (3.18) 


Now using (3.16) and (3.18) we can bound E'{ : 


E E'n 

n>2 ' (fii)" = i \i=l / (ifc)fc fc 

UiSi3a: 

< E E h E (ne-b^-f-'™*) 1 ^ 


e 2 ' |Lfc| 7L fc 


n>2 T tree over (Si)^ =1 \i=l 

{!>•••,«} U iS-Sx 


J-j- e ™(dist hCfii.s,) i) 7Si Sj 

(hj)er 


(3.19) 


where in the sums we keep implicit that (5',)" =1 € SPS^a- 
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Substitute into (3.19) the entropy bound 15 (3.12). Since U,S, 9 x, but not 
necessarily Si 9 x, an extra factor n appears. Moreover observe that |S;| is 
even and > 2 (see the definition (2.1)) and disth (Si,Sj) > 2. Then: 

^ 

n>2 Ttree over (si)i=i,...,n (dij)ij^T \®=1 
Si even >2 dij >2 

( /3 ^-™- i ° g4 ) Si _ /3 A j-j e _^ (di ,_ 1)7d 

*=i / (ij)er 




where jd '■= (1 + o(l)). 

Given n > 2 and <5i,...,d n > 1, the number of trees T over the vertices 
{1,..., n} with given degrees deg-^i) = Si Vi = 1,..., n is exactly 16 

(n — 2)! 

((Si — 1)! • • • (S n — 1)! 

if ^2"_i (Si — 1) = n — 2 and zero otherwise. Furthermore the number of edges 
of T is n — 1. Therefore the bound (3.20) leads to 






(3.21) 


The sum over s gives, as /3 —> oo, 


£<^ ! 


- —t?—5 


s>2 


log 4)s 

^ ((5-1)! 


5>1 


^se-^^^ 51084 " 1 ) 3 = 2e 2 4 10 e-^-^(l + o(l)) . 


s> 2 


(3.22) 


The sum over d gives, as /? —1 00 , 

d>2 

= (l + o(l)) (3.23) 

V d> 2 d>2 / 

= ^ + 0 «) (! + «(!)) = ^ (1 + 0 ( 1 )) 

15 The families of regions (Sj )" =1 such that disth (Si, Sj) = 00 for at least one edge 
(i,j) £ T give zero contribution to the sum, therefore we do not need to worry 
about them. 

16 This is an improvement of the well-known Cayley’s formula. 



22 


771 1 O J 

where we used the fact that 1 — e - ^" = 4 2 (1 -f o(l)) (see lemma 

A.5). Substituting (3.22), (3.23) into (3.21), one obtains 


( \ 71 

e-^-^+^ 0 - (1 + o(l)) ) < (1 + o(l)) . 

/ 

(3. 

Assume /ih — A* v > • Then for /3 sufficiently large (3.24) becomes: 

^^ e -/ 3 (^ h - Mv )+/ 3^±^ 2 e -/ 3 ^± 4 (1 + o(1)) 


(3.24) 


< 


= e -/3 2(Mh-/x v )+/3 ^ (i + 0 (i)) . 

9 


(3.25) 


II. Study of the term £2 • 

The ideas are not far from those already seen for £\. We fix a family of 
regions (Si)" =1 and we assume that there exists ( Lk)k G ^.^(UiS)) such 
that UfcLfc 9 x , otherwise the contribution to £2 is zero. Clearly the line L x G 
, 2 ) 4 ( 1 Ji.Si) that contains x is unique. It is convenient to consider separately 
four cases: 

£2 = £' 2 + £'2 + £'2 + E' 2 " ■ 

In £' 2 we assume n = 0, namely Uj<S) = 0; then L x have to be a maximal 
horizontal line of A . In £'f we assume n = 1, namely there is a unique region 
S and L x may have one endpoint on d° xt S and one on d v A or both on 9® xt S l . 
In £'2 we assume n > 2 and L x has one endpoint on Ui<9® xt S) and one on 
d v A or both on the same <9® xt S). In £'."' we assume n > 2 and L x has one 
endpoint on <9® xt S) and one on d^ xt Sj with i ^ j. 

The case n = 0 is easy to deal with. Indeed, since the unique ( L k ) k G 
such that U kLk 9 x is the singleton ( L x ), 

:= E* = e-^l ia:| 7L* < e-? N (l + o(l)) (3.26) 

(L k ) k ^ 

O k L k 33. 


where N denotes the minimum distance between two different vertical com¬ 
ponents of dA. 

When n > 1, by the lemma 3.3 it holds: 


£’ IK 

U k L k Bx 


\L k \ 


7 L k 


< 


4 £» 1st! 


max 
( L k )k 
U k L k 3x 



\L k \ 


7 L k 


(3.27) 


where it is implicit in the notation that (Lfc)fc G ^^(UiS)) . The factor 7^ 
can take two values (see formula (A.9)), both smaller than 1 for /3 sufficiently 
large (uniformly with respect to Lk), since each line Lk must have at least 
one endpoint on Uj3® xt 5) . 
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Now the case n = 1 is also easy to deal with. Indeed, by bounding the r.h.s. 
of (3.27) by 4' s l e~^ L ^ 7 one obtains: 


E 


// 

2 


se^A 


f)|S| 


E 


n- 


-?l L k \ 




UkLk3x 


< ^ 4 | S | e -m|^| 7 ^ ; 


(3.28) 


then observe that \L X \ > disth(S', x) and use the lemma 3.5 for the entropy: 


V" 


sEE 

s>2 d >1 


32 4 4 S 


7 d 


(3.29) 


i (e ^ 2 , —8 Mh + J 

where 7 d := ( + e p 2 

series in s, d, and use 1 — e _ + 
obtain: 


d )(l + o(l)); finally compute the geometric 
= \e ~(1 + o(l)) (see lemma A.5) to 


< f (4 


10 g-/3(Ph-M, 


’)( 


o-fii 


1 - e ™ 


+ o(l)) (l + o(l)) 


2 25 v / 2 

3 


g-/3 (Ph-Mv)+/3 ' j>1 2 2J 


( 1 + 0 ( 1 )). 


(3.30) 


Assume now n > 2 and that L x has one endpoint on U l d‘^ Kt S l and the 
other one on <9 ext A or both endpoints on the same <9® xt SE By introducing an 
extra factor n we may assume that one endpoint is on d^ xt S Fix a family 
of lines (Afc)| =1 € ^Jzf^UiS)) such that U9 x (namely there is a k 
such that Lfc = L x ). We can introduce the graph G = G((S))j, ( Lk)k ) with 
vertices i € {1,..., n} and edges k € {1,... ,p} : the edge k joins the two 
vertices i,j iff the line Lj. has one endpoint on 9® xt S 1 ,; and the other one on 
d° xt Sj. The graph G is connected, hence G admits at least one spanning 
sub-tree T. Notice that the line L x is not part of this tree. Hence, since each 
factor e _ Tl L fcl j Lk is smaller than 1, 


JJe Lk < e ™ |lx| 7 l* JJ e ™ |ifcl 7 Lk 

fe= 1 fceT 

< e “T dist h (Si,*) ^ ^ JJ’ £ - f (dist h (S i ,S i )-l) g 

(ij)er 
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where 7 s,x : = (^ 2 e +e ^ h 2 dlst h(S,a:)^ (i_|_o(l)) and 7 s,s' := (|e & J + 
e -/3 h 2 (d's t h(S ,5 K 1 )) (l -|- o(l)). Therefore: 


max 

(Lfc)fc, U k L k Bx, 

L x from d* xt Si to d v A 
or from a d* xt S 1 to itself 


H e ’ slLkl 'yL k < 


< e -f dist h (S 1 .*) 7g x max If 

’ T tree over 


g-f (disth^i.Sj)-!) ^ s 


(3.31) 


Now using (3.27) and (3.31) we can bound E/,": 




n\ 

n> 2 (5i)? =1 


E’ n- 


> 2 I I 


lL k 


(L k )fc ? LJfc L k k 
L x from Uid; xt Si to d v yl 
or from a d^ xt Si to itself 


_ n ! 

n>2 Ttree over (Si)? •, \i=l 

“ {I,-,n} 


e ~ ^ dist h (Si ,x) ^ 


I] e-f (disth(Si ’^ ) - 1) 7 SitSj 

(id)er 


(3.32) 


<E E 5 E 


Substitute into (3.32) the entropy bound (3.11): 


^"<E E £ E E E (nf 44 '+’‘ Tl,) 

n> 2 T tree over ( s i )i=l,... ,n d*>l ( d ij'*(i,j)eT \i=l / 

{ 1,... ,n} s ^ even >2 d.-.> 2 


n- 




n 


*2 1 ) 


e 2 


(ij)er 


Idij 

(3.33) 


where 7 d '■= (|e /3J + e ^ h 2 1 ^) (1 + o(l)) and 7 ^ := (^e ^2 + 

(1 + o(l)). Observe that (3.33) is identical to (3.20) up to an 
extra factor ^2 d >l e~^ d ” 7 d,, which equals 


^2 e ™ d *id, = 

d,> 1 


e @2 


1 - e 2 


+ o(l)) (l+o(l)) = (l+o(l)) 
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since 1—e ™ = ^e (l+o(l)). Therefore we assume /Uh — fi v > 

and we exploit the inequality (3.25) to bound the expression (3.33): 


r'" < 


{^je -(l + o(l)) 


= 2 ! 2f +^ e -^2( /ih - #iv )+/3+i+H + 0 (!)) . 

9 


(3.34) 


Finally assume n > 2 and that L x has one endpoint on 9® xt S 1 ,; and one on 
d^Sj with i j . By introducing an extra factor n(n— l)/2 we may assume 
that these endpoints lie on <9® xt S) and 9® xt 52 respectively. Fix a family of 
lines (ifc)fc=i e ^’•^(Ui'S’i) such that UfcL fc 9 x (namely there exists k = k x 
such that Lfc = L x ), then consider the graph G = (?((£))*, (Tfc)fc) with 
vertices i € {1,..., n} and edges k € {1,... ,p} : the edge k joins the two 
vertices i,j if the line L *. has one endpoint on 9® xt S 1 ,; and the other one on 
9® xt 5j . G admits at least one spanning sub-tree T that includes the edge 
k x . Therefore 


Y[e L k < Y[e 2 |ifc| 7 L k < 

k— 1 fceT 

< g-T ( dist h(S , i. a: )+ dist h(S2,a;)-l) ^ ^ TT g-f^disth (Si,S.,)-l) ^ ^ 

(tiler 


where js,S',x ~ (le~P J + e-P^( disth(S ’ x)+disth(s '’ x) - 1 )) (1 + o(l)) and 
IS,S' ■■= (\e~P J + g-/3^( d isth(+S')-i)) (1 + 0 (!)). T h us: 


max 

(Z/fc)/e, Uf~Lk3x 
L x from aj xt 5i to d* xt S 2 



k 


< 


< 


-^(dist h (Si,a;)+dist h (52.a;)-l) ^ 

C Jo 1,02,x 


max 

Ttree over {l,...,n} 

T9( 1,2) 


n 


g-TT (dist h (Si,Sj)-l) 


(tiler 

(id)#(l,2) 


7S, ,5, • 


(3.35) 
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Now using (3.27) and (3.35) we can bound E,"": 

E ('ik ( s “~ t)|s ' iV 


n\ 

n> 2 (5i)^ =1 


E' 


n 


e ” |z " l| 7L fc < 


(^ k ) fc 5 U kL k 3 x 

L x from a d* xt Si to a d* xt Sj with 


s E E 


n ( n - x ) ^ jJJ e -(' 3iai ^-f- lo « 4 )l s ‘l-' 3J 


(5,)" =1 V=i 


n>2 7~tree over {l,...,n} 

T9(l,2) 

. „-T(dist h (Si,s)+dist h (S 2 ,a;)~l) - 


n e ™ (disth(,s ‘’ s,j) i} 7Si,5j 


(i,j)er 

(*7)V( 1.2) 


(3.36) 


Removing the edge (1, 2) from the tree T one obtains two disjoint trees 7i, 7j>. 
By applying to each tree the entropy bound (3.11), one finds: 

#{(5))" =1 G | disth(*Si,a;) = di, dist h (5’ 2 , x) = d 2 , \Si\ = si Vi, 

disth {Si,Sj) = dij V(i, j) G T\ (1,2)} = 

= II H( S i)ier t G 9ST a | dist h (5 t ,a:) = d t , |S;| = s, Vi G 7), 

4=1,2 dist h (Si, S,-) = d^ V(i, j) G 7)} < 



then substitute this entropy bound into (3.36) and obtain: 

n{n — 1) 


< E E 


n>2 7~tree over {l,...,n} 
T9(l,2) 


2 n! 


E E E 

(«*)i=l,...,n d\,d, 2 >l 1,2) 

Si even >2 da>2 


(n f 44s - 4 egr(i) j 

. e -f(di+d 2 -i) 7di+d2 


^ n / 

lb (3 ' 


--^-log4)Si-^J 


(3.37) 


n 


e ™ {diJ 1} 7 d. 


(ij)er 

(hi)V( i,2) 


where 7^ := ( \e~^ J + (d- 1 )^ (1 -p o(l)). 

As already seen, given n > 2 and <5-|,..., <5 n > 1, the number of trees T 
over the vertices {1,..., n} with fixed degrees deg T (i) = Si Vi = 1,..., n is 
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bounded by • Furthermore the number of edges of T different 

from (1,2) is n — 2. Therefore the bound (3.37) leads to: 



The sums over s, d have been already computed in (3.22), (3.23) respectively; 
the sum over di,d 2 gives, as /3 —> oo, 


e -f(<h+*-i ) 7di+d2 

dl,d2>l 


0-0 J 


X 1 — e *Y 


+ 0 ( 1 ) (1 + 0 ( 1 )) 


= 2e /3(Wl+2J) (1 + o(l)) . 


(3.39) 


Substitute (3.22), (3.23), (3.39) into (3.38) and obtain 


Z'l" < 


^ f (l + o(l))) e^ J (1 + o(l)) . (3.40) 

n> 2 V / 


Assume Mh — A*v > IJh + J . Then for /? sufficiently large the (3.40) becomes: 


r"" < 


e ^(l + 0 (l)) 


2 52 £ 


(3.41) 


g—/3 2(/ih —Pv)+^(Mh+2J) ^ o(l)) 


In conclusion, by using the estimates (3.17), (3.25), (3.26), (3.30), (3.34), 


(3.41), and the fact that m = e~ p 5 (l + o(l)) (see lemma A.5), if we 
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assume — p v > j~ J , we find that: 


1 

m 


v 

Pe&A 

supp P9:r 


£(P) e' 


«(p) _ 


= (x; + X" + X' + X" + X'" + X"") (1 + 0(1)) 

< + 2^4 e _ /32(/ih _ Mv)+/3 i i h+2J 


— e 
m 


-^N 


2 25 ' 5 O/ x fl 2/i h +5J 2 52,5 e 4 oo/ Ufl 3Mh+6J 

__ e -/3(Mh-Mv)+/3—^2- _i_ __ e -/5 2(/x h -^ v )+/3-^- 


+ X^ e -/32(Mh-Mv)+/3^pj (1+0(1)) 

= f -e-^ N + (1 + 0(1)) 

\ TO 3 / 


(3.42) 


where IV is the minimum distance between two different vertical components 
of dA and o(l) —> 0 as /3 —» oo (uniformly with respect to N). 

Now we assume that /x v + yy < 0. Thus there exists /?o > 0 such that for 
all /3 > Po the function 1 + o(l) on the r.h.s. of (3.42) is < 2 and the term 
2 __ g/^Ov+^l < 1/32 . There exists 17 also No(/3) such that for all N > No(/3) 
the term X < 1/32 . Therefore if /x v + < 0 (which entails also the 

previous condition /+ — fi v > /th 2 *~ J ), then the inequality (3.42) implies that 

E m* {P) < f 

pe&° A 

supp P3x 


for /3 > f3 0 and N > N 0 (/3). This concludes the proof. 


□ 


4 Proofs of the Liquid Crystal Properties 

In this section we will finally prove that the model behaves like a liquid 
crystal, as stated in the section 1, by means of the cluster expansion results 
obtained in the previous sections. 


4.1 Proof of the theorem 1.2 

We will prove the inequality (1.10) for f\ tX . That one for f x , x can be 
proved analogously; then (1.11) and (1.12) follow since f x = f\ tX + f IjX . 


No = — log — 

u m ° m 
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Observe that 

Z h 

/ £ \ h A\x 

\«/l,x/yl ryYi ? 


where ■Z’ 1 ! is the partition function over the lattice A \ x with horizontal 
boundary conditions including a left-dimer at the site x. Since N > N 0 (f3) 
and disth (x, dA) > N 0 (f3), both partition functions satisfy the hypothesis of 
the corollary 3.2. Hence by the cluster expansion (3.4) the partition functions 
rewrite as 

Z\ = Ca exp ( E* u A ((p t ) t ) 

\p t )tev&>A 

= C A \ X exp ( E u A\x({Pt)t) 

\Pt)t€V3> A \ v 




By applying the definition (2.16), 


CA\x br,x— (1,0) ^l,a:-t-(l,0) 

~c7 = Ai 


Now consider a polymer P £ P A U SZ A \ X . Keeping in mind the definitions of 
polymer (2.12) and polymer activity (2.17), observe that 18 


if disth(suppP, a;) >1 => P £ @> A li &\\ x , g A (P) = g A \ x (P) ■ 
Therefore: 


E* u A\ x {(p t )t) - Y* MWt) > 

> - E* \ u A\*(( p t)t)\ - E* MW 

(Pt)te^A\x (PAtev&A 

Bt: dist h (supp P t , x)<l 3 1\ disth(supp P t ,x)< 1 


And by the inequalities (3.6) and (3.2) applied to both Z\, Z } 1 , 

E* \M( p t)t)\ < Y e( p )e a(p) < 3ji 

(Pt)te^A Pe&A 

Bt: disth(supp Pt, x)<l disth(supp P,x)< 1 


E* |^\x((^)t)| < E Q(P)e a{P) < 2 7 ^ . 

(Pt)t€V&>A\x Pe^A\a, 

3 1: disth(supp P t , x)<l dist h (supp P,x)<l 


18 The condition disth (supp P, x) > 1 guarantees that suppP C A \ x and that 
the polymer P does not include any line Lk having one endpoint oni± (1, 0), nor 
any region Si containing these points. 
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In conclusion one obtains: 

Z A\x ^ K,x-(i,o)bi x+ (i,o) _ ( K rn 

z\ - Ai exp r 5 ¥ 

= \ (l + o(l))) , 

where the last identity follows from the fact that Ai 6 r , a; _(i i o) &i,s+(i,o) = 
El 1] B TtX _ {m B lM0) E^ = ^=(l-§(l+o(l)))^=(l-§ (l + o(l))) (by 
lemma A.8, since there is a left-dimer fixed at x according to Z \\x)' — 

1 + f (1 + o(l)) (proposition A.3), and e -W 8 = 1 - § a & (1 + o(l)) (lemma 
A.5). Finally, since o(l) —> 0 as ft —> oo and o(l) does not depend on the 
choice of x and A, one may obtain the desired inequality eventually increasing 

Po ■ □ 


(/i,*>a = 


4.2 Proof of the corollary 1.4 

Set (fiA,N 0 '■= #{ x e A | disth(x, dA) > N 0 } / \A\. By the theorem 1.2, 
bound (1.11), using also f VtX < 1 — fh,x, one obtains: 

(Z\ orie nt.)A = 4 E - (fv,x)\) > <PA,N oW (l - 4 e"^ ) . 

' ' xeA 


On the other hand: 


Va,n 0 > min <Pl,n 0 

L maximal 
horiz. line of A 


min 

L maximal 
horiz. line of A 


\L\ - 2 Np(p) 

\L\ 


1 _ 9 W) 

N 


□ 


4.3 Proof of the corollary 1.5 

Set ifiA,N 0 '■= #{z € A | disth(x, dA) > N 0 } / \A\. By the theorem 1.2, 
bound (1.12), 

2 _^ ^h_|_ j- 

| (Atransl.) A \ — WjJ E! | {fr,x)A~ A | ^ l fiA,N 0 (/3)^ ^ 2 “+- 1 ¥ 7 ^1,JV 0 ( / S) • 

I I xG2l, 

x ^ even 


On the other hand we have already observed in the proof of the corollary 1.4 
that <pa,n 0 > 1 — 2N 0 /N . □ 
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4.4 Proof of the theorem 1.3 


We will prove the inequality (1.13). (1-14) and (1-15) can be proved anal¬ 
ogously. First of all observe that, since 0 < f\. x . f\. y < 1, 


|(/i,x Av)a - (fiA A (Av)a\ < !og 
Now observe that: 


(fl,x f\,y)\ y (f\,x)\ (Av)\ 

(Ax)\(Av)\ 


(4.1) 


(fl,x f\,y) A ~ 


yh 

^A\x,y 

yh 


(A x)\ = 


yh 

^ A\x 
yh 


(Av)a = 


yh 

h _ Z A\y 
yh 
A A 


where Z\ x , , Ah y are the partition function respectively over the 

lattices yl\x , A\y, A\x,y , with horizontal boundary conditions including 
a left-dimer respectively at the site x , at the site y , at both sites x, y. 
Therefore 


yh yh 
^ A ^ A\x,y 

(Ax) A (Ay) A Z A\x Z A\y 


(Ax Ay)A 


(4.2) 


Since N > N 0 (f3 ), disth(x, dA) > N 0 ((3 ), dist h (y, dA) > N 0 ((3 ), dist h (x, y) > 
No(f3), all four partition functions satisfy the hypothesis of the corollary 3.2. 
Hence by the cluster expansion (3.4) the partition functions rewrites as 


= C A exp ^ E* U A ((P t )t)) , 

\Pt)t€V^>A ' 

Z A\x = C A\X exp f U A\x((Pt)t)) 

VP.l. c^<^03> . , / 


\Pt) t e^A\ x 


Z A\y = Ca\v exp f E U A\y((Pt)t)] 

VD.V ' 


(4.3) 


Z A\x,y CA\x,y e xp 


(Pt)teV&A\y 

E* u a\xA( p a) 

(Pt)teVS>A\x, y 


By applying the definition (2.16), it holds 

Pa PA\ x .y 


PA\x PA\y 


= 1 . 


(4.4) 


Now consider a polymer P G 3? a U &a\x U &A\ y U &A\x,y ■ Keeping in mind 
the definitions of polymer (2.12) and polymer activity (2.17), observe that: 


if disth(supp P. x) > 1 , disth (supp P. y) >1 => 

P G A n &A\x n &A\y n ^A\x,y ) Qa(P) = QA\x(P) = QA\y(P) = 6A\x,y(P) 
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and that 19 : 


if disth(supp P. x) < 1 , disth (supp P, y) >1 => 

P G n & A\y ) \ {&A\x U £*A\x,y) , Qa(P) = QA\y(P) Or 

P G (&A\x n & A\x,y) \ (3*A U ^A\y) , QA\x(P) = QA\x,y{P) Or 

P G £?a n &A\x n &A\y n ^A\x,y > 6a{P) = 0A\y(P) > QA\x{P) = QA\x,y(P) 


and the case disth (supp P, x) > 1, disth(suppP, y) < 1 is clearly symmetric 
to the previous one. Therefore: 


E* - E* u A\x((p t )t) + 

(Pthev&A (Pt)te^ A \x 

E U A\y((Pt)t) + E U A\x,y((Pt)t) < 

(Pt)t&g&A\y (Pt)te^A\ x ,y 

- E \UA({Pt)t)\ + E \UA\x((Pt)t)\ + (4.5) 

(Pt)te^A (Pt)tev&>A\* 

3 1: disth (supp Pt, *)<1 3t: dist h (supp P t , x)<l 

3t': disth (supp P t ' , !/)<l 3t': disth (supp P t i , y)<l 

+ E \ U A\y{(Pt)t) I + E \ U A\x,y((Pt)t)\ ■ 

{Pt)te^ A \y {Pt)teV&>A\x, V 

3t: disth (supp P t , x)<l Eli: disth (supp Pti x )< 1 

3t' : disth ( SU PP Pt' i y) < 1 Eli 7 : disth ( SU PP Pt '» y) < 1 


It is crucial to observe that given a cluster ( P t ) t G ^SPa , since U t suppP t 
have to be connected in Z 2 , 


dist Z 2 (x,y) < dist Z 2 (U t supp P t ,x) + E I SU PP Pt | -1 + dist Z 2 (U t supp P t , y) . 

t 


19 The first possibility, namely P polymer only of the lattices that contain x, 
happens when suppP 3 i or P includes a region S'; containing x — (1,0) . The 
second possibility, namely P polymer only of the lattices that do not contain x, 
happens when P includes a line Lk with one endpoint on x ± (1,0) . The last 
possibility happens when P includes a region S', containing x + (1, 0) (and does not 
verify the other conditions). 
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Hence, assuming that clist^ (U t supp Pf.x) < 1, dist ^2 (U t supp P t . y) < 1, it 
follows 


n« p -) = 


n 


n t \y t \ 


exp — /3 


Mh - Mv 


pt 


Y,\St\-mY,\Lk\^PJ 


n t 


i—l 


k =1 


= exp ( ~ £|suppP*| 


n 


n t \p t \ 


exp — (/3 


j/^ii-Atv m 


-t)Ei s j4En-^ 


Pt 


n t 


i=1 


k=1 


< 


exp (“ (dist Z 2 (a;,y) - i)) 


where P t = (Lfc)|* =1 ) for all t and g*(Pt) is defined as the factor 

appearing in the product over t at the penultimate step. By defining a*(P) := 
suppP|, we have that p*(P) e a *^ is essentially equivalent to g{P) e a ^ : 
we can follow exactly the proof of the theorem 3.1 up to the inequality (3.42) 
and prove that the Kotecky-Preiss conditions (3.2), (3.3) hold also with p* , 
a* and ?n/16 in place of g , a and m/8 (eventually increasing /3o) - Therefore, 
defining C/*((Pt) t ) := u{[P t )t) ^ by the general theory of cluster 

expansion the inequality (3.6) holds also with t/*, g* and a* in place of Ua, 
gA and a. As a consequence: 


E* 1 

(Pt)t&^A 
3t: disth (supp P t , x)<l 
3t / : disth (supp P t /, y)< 1 


< e* K(^)in^) < 

(Pt)t€^A t 

3 1: disth (supp P t , x)<l 
3 1 ': disth (supp P t /, y) < 1 


< e -t(dist z2 ) u ((p t ) t )| n^(p*) 

(Pt)tB^A t 

3 1: disth (supp Pt > :e )<1 
3 1 ': disth (supp P t /, y ) <• 1 


_ x (dist z2 (x,y) — l) E‘ \U*((P t )t)\ (4-6) 

(Pthev&A 

3 1: disth (supp Pt > 

3 1 ': disth (supp P t / , y) < 1 


(3.6) 

< 


e —f(dist z 2 (x,j/)-i) ^ Qi,(P) e “* (p) 

Pe&A 

disth (supp P, x) < 1 


(3.2) 

< e 


(dist z2 (x,y) — 1 ) 3 ™ 

16 ' 
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The same reasoning can be repeated also for the clusters in ^^AXxi 'AS? ' A\y 
and ^^A\x,y Thus, by (4.2), (4.3), (4.4), 4.5, (4.6), one finally obtains: 


(frx fl,y)\ 
(frx)\ (fly )a 


yh yh. 

^A ^A\x,y 


Z \, 


z !l\» 


< exp ^ e -?(dist Z 2 (*,i/)-i) (3 + 2 + 2 + 2) A j . 


The same bound can be shown to hold also for the inverse ratio 
hence by (4.1) we conclude that: 


(fuA liAjdi 

</l,x fl,y)\ ’ 


(fix fry)A 


(fix) A (fry) A 


< g-x( dist z 2(x,y)-l) 


9m 

16 ' 


□ 


4.5 Proof of the corollary 1.6 

Since Aransi. = ttt S xeA - (fr,x ~ frx) > the variance of A rewrites as: 

I yi l xu even 


((Aransl.) 2 )^ ~ ((Aransl.)^)' = ^ 


V. C Xt y 


x h >^h even 


with 

C x ,y ■■= ((/r,x fr,y)\ - (fr ,x)a (fr,y)\) + ((/r,*>A (fry)A ~ (fr,x fry)\) + 

+ ((frx) A (fr,y)A ~ (frx fr,y)\) + ((frx fry)A ~ (frx) A (frvA) • 

By the theorem 1.3, for x, y G A such that disthA dA) > NAB), disth(v, dA) > 
N 0 ((3) and dist h (x,y) > N 0 (P), it holds 

9m 


Hence: 


C xv < 4 — e -f (dist z 2(x,I/)—1) 

x,y - lg 


9 m 


<(Aransl.) 2 )^-((A r ansl.)!\) Z < 4^ £ g -f <dist z2 )+1-^ ^ , 


,y€A 


where we set 

#{(x,y)&AxA | disth(a:, cM) V disth(y, chi') V disth(a;,y)>Afo} 


V’A.A'.ATo : = 


Now observe that 


A \A>\ 


<PA,A,N 0 > min <Pl,L',n 0 > min 

L,Z/ maximal L,Z/ maximal 

horiz. lines of /I horiz. lines of A 


(\L\ — 2No) (\L'\ — ANo) 
\L\\L'\ 


Nn 


> 1 - 2 ^ 1-4^ 


N 


.Nn 


N 
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hence 1 — Ta,a,n 0 < -No /N (6 — 8 N 0 /N ). And on the other hand: 


e ~f ( dist Z 2(*,y)-i) < 


1-11 E 


xGZ 2 

x^O 


(dist z2 (*,0)-l) _ 


= i^iE 4(i( 


-t(d-i) - 


= mi 


d>l 


(1 - e 4 ) 2 


□ 


A Appendix: ID Systems 


Consider a finite line L, that is a finite connected sub-lattice of Z. Consider a 
monomer-dimer model on L given by the following partition function: 


= ^ e-Pn^ e J i( a *i) e -fr(a* r ) 

<«€@I 


denotes the set of monomer-dimer configurations on L (allowing also external 
dimers at the endpoints of L); the Hamiltonian is defined as 

U _ Mh + J „ f sites of I/withl , £_ »(sites of L with dimer but! 

L 2 monomer J 2 * l neighbor to monomer in L J 

xi,x r denote the left and the right endpoint of L respectively; I\,I T represent the 
interaction among the configuration on L and the boundary condition outside its 
endpoints. 

This one-dimensional system can be described by a transfer matrix T over the 
three possible states of a site, l = “left-dimer”, r = “right-dimer”, m = “monomer”: 

/ T(l,l ) T(l,r) T(l, m) \ (0 1 y/aF\ 

T = T{r,l) T{r,r) T{r,m) ] := 1 0 0 , (A.l) 

\T(m,l) T(m,r) T(m,m)J \0 y/ab a J 


o + J 

where to shorten the notation we set ya := e p 3 the transfer contribution 

of a monomer 20 , yfb := the transfer contribution of a site with a dimer 

but neighbor to a monomer. Two vectors are also needed to encode the boundary 
conditions: 


B\ = (Bi(l) Bi(r) Bi(m)) := {e™ e /lW y/^e 1 ^) , 


( Br(l) \ 

( \ 

(A.2) 

B r = B T (r) 

:= e ,A - r) 


\B r (m)J 

\yfie 1 '™) 



Proposition A.l The partition function of the system rewrites as a bilinear form: 

Z L = B\ T |i|_1 B r . (A.3) 


20 The transfer energy of a monomer is half the energy of a monomer because it 
appears during two “transfers”. 
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Proof According to the previous definitions it is clear that for every configuration 
a £ {l, r , m}^ 

t(a £ 3> l ) e~ pHL{cx) = 

= ^L 1(ai=rn) T(ai,a 2 )T(a 2 ,a 3 ) ... T(a lLpi ,a lLl ) . 

Therefore 

Zl = ^2 Bi(ai) T(ai, 02 ) T(oi 2 , 03 ) • • • T(a|n|-i, cx\l\) B r (ot\ L \) 

at£{l,r,m} 1^1 

= B\ T 11,1-1 B r . 

□ 

Assume for the moment that the transfer matrix T is diagonalizable. Denote 
by Ai, A 2 , A 3 its eigenvalues and by E ^, E, Ej 3 ' 1 , e [ B , E j (2) , E[ 3 ^ the correspond¬ 
ing right (column) eigenvectors and left (row) eigenvectors, normalized so that 
E^E^ = 1 for i = 1,2,3. 

Corollary A.2 

Z L = J2 B i E x ] Br ■ (A.4) 

*=1,2,3 

Proof Since we are assuming that T is diagonalizable, it holds T = P D P^ 1 
where D is the diagonal matrix of eigenvalues, P is the matrix with the right 

eigenvectors on the columns, P~ 3 has the left eigenvectors on the rows. Then 
= pd \l\-i p-i and 

3 

B\ T |l|_1 B t = (B\ P) D |l|_1 (P _1 B r ) = ^(Bi£ r (i) )A [ i|-1 (E^B r ) . 

i=1 

□ 


Now our purpose is to diagonalise the transfer matrix T when (3 is large. 


Proposition A.3 For all j3 > 0 the transfer matrix T is diagonalizable over R. 
Its eigenvalues are 


Al = 1 + y (1 + 0(1)) 


A2 = -1+y (1+0(1)) ( A - 5 ) 

A 3 = a — ab — a 3 b (1 + o(l)) 


as (3 — > 00 . 

Proof The eigenvalues Ai, A 2 , A 3 are the (complex) roots of the characteristic poly¬ 
nomial of T, that is 

p( A) := det(A/ — T) = -ab + (A - a)(A 2 - 1) . 

For all j3 > 0 it turns out that p has 3 distinct real roots 21 , hence T is diagonalizable 
over the reals. 

21 The discriminant of the cubic is A = 18a(l—&) + 4a 2 (l— 6 ) + a 2 + 4 — 27a 2 (l— b), 
which is strictly positive for all 0 < a, b < 1 , (a, b ) ^ ( 1 , 0 ). 
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As /3 —> oo, p( A) —> A(A 2 — 1) hence Ai —> 1, A 2 —> — 1, A 3 —> 0. Thus it is 
convenient to write Ai = 1 + £1 , A 2 = — 1 + £2 , A 3 = a + £3 with £, —> 0 as —> 00 
for i = 1,2, 3. Now expand the polynomial p in powers of £i and truncate it at the 
first order: 

0 = p(Ai) ~ —ab + (1 — u + £ 1 ) (2ei + £ 1 ) = —ab + 2£i (1 + o(l)) 

ab ,,,, 

=> £1 = y (! + °(1)) ; 

0 = p( A 2 ) ~ —ab + (—1 — a + e) (— 2£2 + £ 2 ) = —ab -f- 2£2 (1 + o(l)) 

ab ,,,, 

=> £ 2 = y (1 + 0 ( 1 )) ; 

0 = p(As) = -ab + £3 ((a + £3) 2 — l) = —ab — £3 (1 + o(l)) 

=> £3 = —ab (1 + o(l)) . 

In order to find the following order of A3, now one can write A3 = a — ab (1 + £3) 
with £3 — > 0 as (3 —> 00 and repeat the procedure: 

0 = ^ = 1 + (1 + £3) ^ a 2 (1 + o(l)) — 1 ^ = a 2 (1 + o(l)) —£ 3(1 + o(l)) 

=> £3 — a 2 (1 + o(l)) • 


□ 


Proposition A.4 The right eigenvectors of the transfer matrix T are 

1 Z 1 “ I (1 + °( 1 ))\ 

E r (1) = 4= 1-f (l + o(l)) 

v 2 y Vab (1 + o(l)) J 

1 / 1 + f ( 1 + °( 1 )) \ 

S r (2) = — -1 - | (1 + o(l)) (A.6) 

\ Vab (l + o(l)) J 

(— a\fab(l + o(l))\ 

E™ = —Vab (1 + o(l)) 

V 1 + a (1 + o(l)) / 


and moreover 

E ( r 2 ) (l) + E?\ 2 ) + V^bE?\ 3 ) = ^ (l + o(l)) 
£ r (3) (2) + VabEi 3 ) { 3 ) = —a 2 Vab (1 + o(l)) 


as /3 —y 00 . The left eigenvectors are obtained by a simple transformation: E = 
a(Er^) for i = 1,2,3, where 


( y i\ 

a I V2 I := (V2 Vi V3) . 

\V3 J 

Proof The right eigenvectors E r associated to the eigenvalue A are the non-zero 
solutions of the linear system 

ZA(A - a)\ 

(XI - T) Ey = 0 E r = A -a U, (el. 

\ Vab J 
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And the left eigenvectors E\ associated to the same eigenvalue A are the non-zero 
solutions of the linear system 

E\ (A I — T) = 0 +> Ei = (A — a A(A — a) \/ab) t , ( £ K . 

The desired normalization E\ E r = 1 can be obtained by choosing 

t = W2A(A — o) + ab 

in both cases. Now to conclude the proof it is sufficient to exploit the estimates of 
the eigenvalues given by the proposition A.3. □ 

The formula (A.4) together with the estimates of propositions A.3, A.4 give us 
a complete control on the one-dimensional system on L at low temperature, for 
every choice of the boundary conditions. 

We concentrate on providing an estimation of the quantity Rl defined by (2.9), 
since it is needed in the section 2. We have to distinguish three cases, according to 
where the endpoints of L lie. 

transfer matrix T are 
—a + ab (1 + o(l)) 

| it holds 

as /3 —> oo . (A.7) 

A.3. □ 

Lemma A.6 If xi £ <9® xt Sj, then as 0 —> oo 

B,£ r (1) = ^ (l + o(l)) 

B,S r (2) = -^ (l + o(l)) 

Bi£ r (3) = v'o (l + o(l)) • 


Lemma A.5 The ratios of the eigenvalues of the 

^ = -1 + 0601+0(1)) , ^ = 

as fl —> oo. In particular setting m := — log | A 2 /A 1 

e" m = 1 e 3 •' (l + o(l)) 

Proof ft follows immediately from the proposition 


If x r £ dr'Sj, then the same estimates hold for Ef- 1> B r , E[ 2 ^B v , E^B T respec¬ 
tively. 

Proof If x\ £ <9® xt Sj then by (2.7) and (A.2) the vector describing the boundary 
condition on the left side of the line L is B\ = (O Vh v / a ) . Then the estimates for 
B\ i = l,2,3, are computed using the proposition A.4. □ 


Lemma A.7 If x 1 £ <9i A, then as /3 —> 00 


B\ £ r (1) = 


B\ E, 


f ^75 (l — f (1 + o(l))) if the h-dimer on x\ — ( 1 , 0 ) has fixed position 
] (l — | (1 + o(l))) if the h-dimer on x\ — (l, 0 ) has free position 

— (l + | (1 + o(l))) if the h-dimer on x\ — ( 1 , 0 ) is fixed to the left 


B\ E^ = i (l + f (1 + o(l))) if the h-dimer on xi — (1, 0) is fixed to the right 

(1 + o(l)) if the h-dimer on x\ — ( 1 , 0 ) has free position 

—a 2 \fab (1 + o(l)) if the h-dimer on x\ — ( 1 , 0 ) is fixed to the left 
—a\/ab (1 + o(l)) if the h-dimer on x 1 —( 1 , 0 ) is fixed to the right or free 


( 3 ) _ 


If x r £ d r A, then the same estimates hold respectively for e[ 1] B v , E[ 1 ' 1 B r , E\'^ B t 
after substituting: x\ — ( 1 , 0 ) by x r + ( 1 , 0 ), “left" by “right” and “right” by “left”. 
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Proof If Xi £ d\A then by ( 2 . 7 ) and (A. 2 ) the vector describing the boundary 
condition on the left side of the line L is: B\ = (0 1 Vab) if a left-dimer is fixed 
on *i —( 1 , 0 ); B\ — (1 0 0 ) if a right-dimer is fixed on xi —( 1 , 0 ); B\ = (l 1 i/ab) if 

on xi —( 1 , 0 ) there is a h-dimer with free position. Then the estimates for B\ E^\ 
i = 1 , 2 , 3 , are computed using the proposition A. 4 . □ 


Lemma A.8 If x i £ diA, then as /3 —> oo 


j 775 (1 — f (1 + °(1))) */ the h-dimer on x\ — (l, 0) has fixed position 
1 \[2 (l — | (1 + o(l))) if the h-dimer on x\ — (1, 0) has free position 


{ — (1 + | (1 + o(l))) if the h-dimer on x\ — (1, 0) is fixed to the left 

(l + § (1 + o(l))) if the h-dimer on x\ — (1,0) is fixed to the right 

^7= (1 + o(l)) if the h-dimer on xi —(1, 0) has free position 

(3) I — a 2 \/ab (1 + o(l)) if the h-dimer on x\ — (1, 0) is fixed to the left 
1 11 = < — 

I —ayab{ 1 + o(l)) if the h-dimer on xi —(1,0) is fixed to the right or free 

If x r £ d r A, then the same estimates hold respectively for B v , E [ 2 ' 1 B r , £)+' B T 
after substituting: x\ — (1,0) by x r + (1,0), “left" by “right” and “right” by “left”. 


Proof If Xi £ d\A then by ( 2 . 7 ) and (A. 2 ) the vector describing the boundary 
condition on the left side of the line L is: B\ = (0 1 \fab ) if a left-dimer is fixed 
on xi —( 1 , 0 ); T?! = (10 0 ) if a right-dimer is fixed on xi —( 1 , 0 ); B\ = (l 1 \faS) if 

on xi —( 1 , 0 ) there is a h-dimer with free position. Then the estimates for B\ Ei l \ 
i = 1 , 2 , 3 , are computed using the proposition A. 4 . □ 


Proposition A.9 Denote by o(l) any function u>(/ 3 , /Xh, J) that goes to zero as 
ft —> oo and does not depend on the choice of the line L nor on A. Then for every 
line L £ Jzf/t(U(«S'j), Sj £ S?a pairwise disconnected, A CZ 2 finite, it holds 

\Rl\ < (A.8) 


where the quantity 7 l can be chosen as follows: 

( + e- /3 ^ |i| ) (1 + 0 ( 1 )) if XI £ U id^Si , Xr £ Uidr*Si 

7 L ■■= < 


^(1 + 0(1)) 

ifx 1 £ Uid r ext S;, x v £ VJidyA 
or vice versa xi £ U id\A , x r £ U 

,1 + 0(1) 

if xi £ diA , x r £ d r A 

(A. 9 ) 


Proof • Suppose Xi £ d^ xt Si and x r £ df xt Sj . The definition ( 2 . 9 ) and the corollary 
A .2 give 


XiRl = -rrfy - BxE^E^B, 

= ) BiE[ 2 ) e[ 2 ) B t + BiE^ e[ 3) B Y . 

\Ai / \Ai J 
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By the lemma A.5 | A 3 /A 1 1 < a | A 2 /A 1 1 when /3 is sufficiently large. Therefore, using 
also the estimates of lemma A. 6 , one finds 


A2 

Ai 

+1-1 ^ 

A2 

Ai 

+1-1 ^ 


+ a |L ) (1 + o(l)) . 


\Rl\ < 


• Suppose now *1 £ 9® xt 5j and x r £ d r A. The definition (2.9) and the corollary 
A.2 give 

Z L 


Ai /2 7?l = 


x \l\- i e ( 1) Bi 


- BiE£ 


(i) 


A 2 V L|_1 BiE^E^B, + {\ 3 \ L '- 1 BiE^ e[ 3 ) B? 


Ai 


E[ 1] B r 


Ai 


e[ 3) b t 


By the lemma A.5 | A 3 /Ai | < a | A 2 /Ai| when /3 is sufficiently large. Therefore, using 
also the estimates of lemmas A. 6 , A. 8 , one obtains 


1^1 < 


\L 1-1 


with 


if fixed h-dimer on a; r + (l,0) (1 + o(l)) + a) L ^ 1 0(a 2 Vb) 

if free h-dimer on x r + (l, 0) ^ (1 + o(l)) + a W- 1 (1 + o(l)) 

f ^(l + o(l)) ,/h 

d+od)) ^ (i+o(i)) - 

• Suppose now x\ £ d\ A and x T £ d r A. The definition (2.9) and the corollary A.2 
give 

Z L 


Rl = 


- 1 


A 2 y i|_1 BiE^E^B, + f X 3 \ L '- % BiE^ 3) E^ 3) B r 


A ^^BiE^E^B 

A 2 

Ai 


BiE[ 1] E[ 1] B, 


Ai 


BiE^E^B, 


By the lemma A.5 | A 3 /Ai | < a |A 2 /Ai| when /3 is sufficiently large. Therefore, using 
also the estimates of lemma A. 8 , one obtains 


\Rl I < 


\L 1-1 


with 


7 = S 


if fixed h-b.c. on both sides 1 + 2 a (1 + o(l)) + a) L 1 0{a 3 b ) 

if fixed h-b.c. on one side, vr (1 + o(l)) + a ^ _1 0{a 3 b) 
free h-b.c. on the other one 4 

[if free h-b.c. on both sides (1 + o(l)) + a '^ -1 ^ (1 + o(l)) 

f 1 -{- 2 a (1 + o(l)) 

f (l + o(l)) < l + o(l). 

^ + »SP) (1 + 0(1)) 
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B Appendix: Cluster Expansion 


In this Appendix we state the main results about the general theory of cluster 
expansion used in this paper. The condition that we adopt to guarantee the con¬ 
vergence of the expansion is due to Kotecky-Preiss [11]. For a modern proof we 
refer to [17]. 

Let ^ be a finite set, called the set of polymers. Let g : 3? —> C, called the 
polymer activity, and 8 : 3 s x 2? —> { 0 , 1 }, called the polymer hard-core interaction, 
such that 8(P,P) = 0 and 8(P,P') = 8(P',P ) for all P,P' £ £?. Consider the 
polymer partition function: 

z ■■= e n ^ n 5 ( pp/ ) 

pes?' pp'eg?' 

P ’ (B.l) 

= zl e iW«) n *(*>*) ■ 

q >0 Pl,...,P q e3 /> t= 1 t<s 


A family of polymers (Pi,... ,P q ) is called compatible if 5(Pt, P s ) = 1 for all t s ; 
otherwise it is called incompatible. Observe that in the partition function Z only 
the compatible families of polymers give non-zero contributions. 

A family of polymers (Pi,... ,P q ) is called a cluster if the graph with vertex set 
{ 1 ,..., q} and edge set {(t, s) \ 5(Pt, P s ) = 0} is connected. 

Theorem B.l Suppose that there exists a: & —> [ 0 , oo[, called size function, such 
that the Kotecky-Preiss condition is satisfied, namely: 


Y \Q{P)\e a(P) < a(P*) V P* £ . 

pe& 

S(P,P *)=o 


Then: 


where the 


i°g2 = e i e 

<?>o p 1 ,...,p q e&> \t =i / 

series on the r.h.s. is absolutely convergent and 


u(Pi,...,P q ) := Y. (-1) |B| - 

G=(V,E) connected graph 
V={l,...,q} 

EC{(t,s)\S(P t ,P s )=0} 


(B.2) 


(B.3) 


(B.4) 


Moreover, for all £ C 3? 


E 

q >0 H P 1 ,...,P q e^ > 
3 t : Pt £ $ 


n«(*: 


\u(Pi,. 


-,P ,)I < E \e(.P)\z a(P) ■ 

P£&> 

peg 


(B.5) 


It is worth to observe that if (Pi,..., P q ) is not a cluster then u(Pi ,. .. , P q ) = 0. 
Therefore only the clusters of polymers (that are infinitely many) give non-zero 
contributions to the expansion (B.3) of logiL 
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